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Abstract 

We show that an escape from the potential minimum of Fabry-Perot interferometers can be detected 
measuring the associated sudden change of reflectivity. We demonstrate that the loss of information that 
occurs retaining only the sequence of escapes , rather than the full trajectory, can be very mild and can lead 
to an effective method to reveal the noise intensity or the presence of a coherent signal. 
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I. INTRODUCTION 



Many noisy systems, from condensed matter to nuclear physics, are characterized by a sudden 
passage from an initial state to the final product, the fluctuation-induced escape problem treated 
by Kramer in a celebrated paper [1]. In such cases one measures the escape times (ETs) to retrieve 
information about the physical systems, i.e. the activation energy. To restrict the available data to 
escape events rather than to consider the full trajectory might be a simplification of the measure- 
ment, but it entails a loss of information, in as much as a continuous variable is substituted by a 
discrete sequence of events. While in many systems the full trajectory is practically unavailable, 
i.e. in chemical reactions, in some cases one should balance the advantages of the ETs measure- 
ments with the loss of information. This is the case, as we will show, for a Fabry-Perot (FP) 
cavity. We demonstrate how it is possible to simplify the FP measurement detecting the change in 
reflectivity, at the cost to only retrieve the ETs from the metastable minimum. We also show that 
the loss of information can be lessened by exploiting the whole ETs distribution with a likelihood 
ratio test. 

The letter is organized as follows: first we illustrate the possibility to measure the ETs from a 
minimum of the FP cavity. Then we explicate the underlying mathematical model that describes 
the escape process and that will be employed to generate the ETs distributions to derive the prop- 
erties of the estimates for the noise intensity and the presence of a coherent sinusoidal are outlined. 
In the final section we conclude and discuss. 

We demonstrate how it is possible to simplify the FP measurement detecting the change in 
reflectivity, at the cost to only retrieve the ETs from the metastable minimum. We also show that 
the loss of information can be lessened by exploiting the whole ETs distribution with a likelihood 
ratio test. 

The letter is organized as follows: first we illustrate the possibility to measure the ETs from a 
minimum of the FP cavity. Then we explicate the underlying mathematical model that describes 
the escape process and that will be employed to generate the ETs distributions to derive the prop- 
erties of the estimates for the noise intensity and the presence of a coherent sinusoidal are outlined. 
In the final section we conclude and discuss. 
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II. FABRY-PEROT CAVITY OSCILLATIONS: PROBLEM STATEMENT 



A FP cavity is formed by suspended mirrors that make up an optical cavity for a high power 
laser. The pendulum provides a mechanical restoring force, while the spatial period of the laser 
radiation pressure induces multistability [|2|-|4|]. The FP constitutes the latest technology in pre- 
cision measurements of weak signals Q5D also in the quantum regime [6D or gravitational waves 
detection [7]. A weak perturbation forces oscillations of the mirrors that can be detected by the 
change in the laser power accumulated in the optical cavity. However, noise (included the noise 
that enters the system through the measurement devices) should be constrained to reveal the very 
weak interaction of gravitational waves with matter, for noise can induce spurious oscillations of 
the pendulum that mask the effect of the gravitational waves. The key feature that we want to 
exploit is the fact that when the displacement of the pendulums matches a semi-integer number 
of the laser light wavelength, a sudden change of reflectivity occurs, see the upper part of Fig. [T] 
It is thus possible to measure the passage across a point close to the maximum of the potential, 
denoted by V + and V\_ in the lower part of Fig. [TJ and to record the time spent in the potential 
minimum. The role of the noise is, in this setup, twofold (as in stochastic resonance [8]): on one 
hand it helps to overcome the potential barrier and on the other hand the noise covers the signal. 
We propose to characterize the FP through the analysis of the time to escape from a minimum 
of the optomechanical potential [9] to exploit the statistics of occasional large excursions (made 
possible by the noise) that can be detected as escape events. The advantage of using Escape Times 
(ETs) stems from the simplicity of the measurement, because the change of reflectivity can be 
detected with the accuracy of optical measurements. It is thus possible to record the time that 
has elapsed between two changes, or to measure an ET from a minimum, with an analysis of the 
light transmitted through the mirrors without mechanical interaction with the pendulums. In other 
words, the measurement of the reflectance allows us to analyze the dynamics of the pendulums 
without the application of the servo system, i.e. to follow the almost unconstrained motion. This is 
a simplification of the measurement technique that could also pave the way for an analysis of the 
quantum regime [|10j]. The simplification has a cost: the ETs abridge the trajectory and therefore 
some information is lost. However, as has already been pointed out in the case of superconduct- 
ing Josephson detectors, the loss is relatively mild and can be partially compensated by a careful 
analysis of the ETs 111 in . 

In fact ETs have proved to be effective in evidentiating subtle effects such as the quantum tun- 



3 



neling from different attractors [U2I1 . the granular nature of the charge fl 1 3M and the discrimination 
between classical and quantum activation II 1411 i n Josephson junctions, or the critical exponents of 
bifurcations in micromechanical oscillators [|15ll . 

Finally, let us remark that if the system is reset after the crossing of the reflectivity minimum, 
the measured time corresponds to a first passage time, while it amounts to an escape if the pen- 
dulums are left to freely oscillate, for the inclusion of an additional backward probability current 
Jl^ll . However, when the absorbing barrier is sufficiently close to the maximum of the potential 
the ET and the first passage time are comparable, and there exist formula for the passage from 
one quantity to the other QJ/ZL ll8|] . In the following we use the escape time and first passage time 
interchangeably, neglecting the difference. 



HI. MODEL EQUATIONS 

To make these ideas quantitative for FP cavities, let us start with the analysis of the steady state 
fields at the mirror points (see the upper part of Fig. [T). The FP cavity stores a circulating power 
P • 

1 arc- 

Pcirc=\H 2 = X 2 ^nA(<P). (1) 

Here x — Phi/Pin = t\/{l — rir 2 ) 2 denotes the FP gain (r 4 and U are the reflection and trans- 
mission coefficients, P M and P in are the maximum and input laser power, A is the Airy function 
A = 1/(1 + J-"sin 2 0), T = 7r v /rir 2 /(l — rir 2 ) the Finesse of the cavity, and the phase of 
the light that also determines the half- width 5(p of the resonance, 5(p = ix j2T . The maximum 
stored power corresponds to the peaks of the Airy function (0 = nir). The detuning phase can 
be changed to match the resonance condition by varying either the length of the cavity L or the 
frequency of the input light. We introduce a further simplification by assuming that the mirrors 
move slow enough that the laser transients can be neglected 

This corresponds to the adiabatic approximation where the electromagnetic transients can 
be neglected, for the distance is large compared to the characteristic time: 2tt/uj >> 
(2L/c) \ ln(rir 2 )| _1 (c is the speed of light and cu is the resonant angular velocity of the pen- 
dulum). In the adiabatic limit the laser builds up the energy in the cavity much faster than the 
mirrors displacements, and the energy stored in the cavity is a function of the positions of the mir- 
rors and not of the previous history. Consequently, the fields inside the cavity are only determined 
by the positions x\ and x 2 of the mirrors through the relative coordinate z = %2 — x\ and the phase 
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FIG. 1: The upper part shows the schematic of the mirrors. The incident electromagnetic 
field ipin is partially reflected (VVe/) and partially transmitted The fields at the mirrors con- 
sequently read Ipi = ti1p in - nV>4, *Ptrans = H^\Z~ % ^, *p3 = ~tp2, Aef = «"l^in + ^1^4- 

In the lower part is shown the potential of the pendular FP interferometer (solid line) as 
a function of the displacement between the mirrors. The dashed line is the reflectivity, to 
be monitored to measure ETs. The minimum of the potential V~ is marked by a black dot, 
while the minima of the reflectivity (indicating the escape) are marked by white dots at the 
potential V+ and V_[. The parameters of the system are: F = 50, = 5, t\ = t% = 0.97, 
C = 10~ 6 . C is the loss factor. 
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of the light cf) = 2n(L + z)/\ (A is the wavelength of the light). The normalized displacement 
z = zj A is governed by the following normalized stochastic differential equation (for details see 

Hi): 

z + 7i = -z + IL m A(2ttz) + e sin(fit + fa) + f (t), (2) 

where overdot denotes differentiation with respect to the time normalized to lOq 1 , ITm = 
Rm/^XcUq is related to the maximum radiation pressure Rm = (2/c)Pm and fi is the reduced 
mass of the mirrors. The friction constant 7 = 7/wo is given by the pendulum dissipation constant 
7 divided by ojq. The parameter that tunes the nonlinearity is the Finesse, that for interferometers 
typically lies in the interval T ~ 1 -f- 1000. In normalized units the pendulums are subject to the 
potential shown in the lower part of Fig. [Q 

v2 



z- tan- 1 [ v / ^TItan(27r^)] +k[2z + 1/2\ 

Viz) — 11m . 

K ' 2 2nVTTT 



(3) 



In the model Eq. © a sinusoidal forcing term of normalized amplitude e and a stochastic noise 
term £ are included. In the upper part of Fig. [T] is also displayed an easily measurable quantity 
introduced and discussed in the previous section: the FP reflectivity (see [7J] for details). The 
reflectivity is a function of the normalized mirrors displacement through the phase of the light cf>: 

t j r2e -2i<t> 



r c (z) :-- 



i>ref 









1 — rir 2 e 



(4) 



The dynamical change of the mirror displacement, determined by the effect of both noise and 
the signal, see Eq.©, affects the reflectivity that abruptly changes in the points indicated in Fig. 
[Q as discussed in the previous section. 

Let us finally discuss an approximation regarding the noise term £(t) introduced in ©. In 
actual pendular FP interferometers the noise is colored, being dominated at low frequency by 
the seismic noise. The frequency spectrum is affected by several other sources of noise such 
as substrate Brownian noise, thermal fluctuations of the expansion coefficient of the test mass, 
and the noise arising from the mechanical dissipation of the mirrors coating [|20ll . However, to 
maintain the discussion at a general level, we assume the total external noise term £(£) Gaussian 
and white: < £(£)£(£') >= 2D§(t — t'). In this setup, in the next sections we focus on two tasks: 
the estimation of noise intensity D and the detection of a sinusoidal signal. 
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IV. NOISE INTENSITY ESTIMATION 



In this section we propose an estimation procedure for noise level based on escape time statis- 
tics. We note that in the absence of an applied signal (e = 0) the system leaves the potential 
minima of Eq.© at the Arrhenius rate oc exp (— jAV/D), for the external noise does not follow 
the fluctuation dissipation theorem 112 ill (AV = V + — VL, where V+ and V- denote the exit point 
and the minimum of the potential, respectively, as shown in the lower part of Fig. [TJ). The pendular 
interferometer is characterized by a very low damping term 7 ~ 10~ 6 Il22ll . thus the system equa- 
tion © should be handled in the extremely low dissipation limit. Some analytical approaches have 
been proposed to extend the vanishingly damping treatment of stochastic differential equations up 



to finite values Q23 



2411 . We have employed the method described by Ref. I125ll to find the ET of 



Eq.© in the absence of the signal, inserting the potential Eq.© in the energy diffusion limit of the 
escape rate I126ll . Using the derivative of the action I in the minimum of potential and integrating 
by parts we obtain: 



< t > 



^ I I(E)e- E ^ D dE 

AV 
~D~ 



+ 



V+ e -E'~i/D 



dE' 



(5) 



that is more convenient for both numerical and analytical treatment. In fact, exploiting the approx- 
imation : 

1(E) ~ (V - V^)2n/u r + (V - 1/_) 2 /(V_)/2, (6) 

where u r is the well resonance, the average ET, Eq.© can be analytically evaluated : 

1 

1 

7pj v\ p, 



< T >-- 



1 + 2 P )<t h > --(1 - 2p) log(l + P1 AV/D) + 

1 



+ -e-V p (l + 2p) 

7 



Ei 



- Ei -fAV/D 



(7) 



Above, Ei(-) is the exponential integral function (see [|27ll for definition and details). For p = 
oj r I(VJ)D / {2,11^) = 0, Eq.© reduces to the ET of the harmonic oscillator : 

< t h >= (Ei( 7 AV/D) - log^AV/D) - je) 1 7, (8) 

where 7^ is the Euler-Mascheroni constant. From the numerical distributions of the ETs shown 
in Figf2]it is evident that ETs distribution can be exploited to determine the intensity of the noise. 
One can compare the average ET with the predictions of Eq.© to estimate the best fit temperature. 
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This amounts to inferring the noise level from the distributions of Fig. [2] with statistical estimation 
of the parameter D in Eq.©. To evaluate the performances of such procedure we have numerically 
retrieved the ETs distributions employing a quasi-symplectic modified velocity Verlet algorithm 
with velocity randomization for the integration of stochastic differential equations [Q]. Moreover, 
we have found consistent results with another leapfrog algorithm [|29ll that has proved to be very 
efficient even at extremely low dissipation BOIl . Finally, the algorithms have been tested against 
the known estimates for the washboard potential at very low dissipation (7 << 1) and noise 
(D « 1) I125ll . We have found, for the whole dissipation range 1CT 6 < 7 < 1CT 2 , that the theory 
lies within the 95% confidence limit of the numerical simulations and that the distortion of the 
temperature estimate is below 5%. 

The efficacy of the method for the estimate of the noise level has been evaluated computing the 
estimator variance a 2 N as a function of the sample size N and of the noise intensity D, under the 
hypothesis that ETs are described by an exponential distribution. This estimator is the Maximum 
Likelihood estimator, and its large N behavior can be determined according to 

^ T »]' 

" ~ NJ(< t >) ' m 

where the function T>(< t >) is the inverse of Eq.© w.r.t. the noise level D, and J(< r >) is 
the Fisher information , that for the exponential distribution reduces to J(< r >) = l/(< r >) 2 
BHD. In Fig. [3] we show that the relative error increases with noise level up to a saturation point 
(around D ~ -yAV) that depends, as expected, upon N^ 1 ^ 2 and exhibits a weak dependence upon 
p (see the inset). The estimate of the noise level improves by increasing 7AV i.e. the dissipation 
and/or the potential barrier. From the physical consideration that a minimum time occurs before 
escape 113211 (see also Fig. |2j), it is evident that the ET density departs from the exponential model 
for short escapes. Such deviations from the exponential distribution lead to an overestimate of the 
relative error (as evident from the fitted curve in Fig. |3), since the numerical result is about 30% 
smaller. 



V. SINUSOIDAL SIGNAL DETECTION 



We shall now focus on sinusoidal perturbations, i.e. e ^ in Eq.©. A sinusoidal signal 
modifies the escape process Il32l - l34ll (see Fig© and could therefore be revealed with an analysis 
of the ET distributions. The analysis can be performed with a direct Sample Mean (SM) procedure 
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Signal (dashed) 
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D = 0.025 
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FIG. 2: Probability distribution of the ET with (dashed line, £ / 0) or with- 
out (solid line, e = 0) an applied signal for two different values of the noise 
intensity (D = 0.025 and D = 0.0125). The signal amplitude reads e = 0.05, 
the phase $o = 0, and the frequency U = 2.6. The parameters of the system 
are: T = 2, U M = 2.1, AV = 0.07, u r = 2.9, 7 = 10~ 6 . 

that considers the average ET. The SM detection of the signal is performed observing a change of 
the average (< r E > for e ^ instead of < r >), a method that has been used to investigate the 



36|]. We underline the role of the noise: 



effect of a periodic signal in Josephson junctions [1351. 
was it be absent, the system would never experience an escape over the barrier, (< r >— > 00 in 
the zero noise limit). Therefore the signal induced displacements are much smaller than the light 
wavelength, as required to make the detection method interesting. 

The SM escape can be optimized with an appropriated choice of the detector characteristics that 
depends upon the noise intensity, and the signal properties to achieve the matching between the 
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FIG. 3: Relative accuracy of the detected noise level as a function of the external noise 
intensity for different values of the sample numerosity. Symbols refer to simulations 
of the Langevin Eq.© for N = 1000, T = 20, U M = 2, AV = 0.042, D = 
0.00042 (filled points) and D = 0.00084 (empty circles). The dashed line denotes the 
fitting. The inset shows the asymptotic behaviors as a function of the parameter p that 
represents the deviations of the model Eq.(f2) from the harmonic oscillator (p = 0). 



noise induced ET and the signal period (see the lower part of Fig. H]), around f2 ~ 10 x ) B4l 13511: 



(10) 



Eq. (TTOl ) ensures a resonant activation, i.e. the possibility to exploit noise to facilitate signal 
detection [8]. To improve the detection, we have employed the Likelihood Ratio (LR) technique 
that analyzes the data exploiting the full distribution of the ETs. LR amounts to performing a 
statistical test to distinguish if the finite sample collected by the experiments is better fitted by the 
unperturbed escape (e = 0) or by an oscillating distribution (e ^ 0). It consists in averaging the 
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samples of a random variable A obtained from escapes by the nonlinear mapping 

A = In [4^1, (ID 
J( T ) J 

where the probability density functions / e (r) and /(r) of the escapes are numerically computed 
with an approach similar to the analysis of the ETs of Josephson junctions y. The LR leads 
to a remarkable improvement of the performances with respect to the analysis of the average ET, 
because the averages < A > and < A e > (computed for e = and respectively) carry 
much more information. This is shown in Fig J4] through the index dxc 0711 . a rough measure of 
the efficiency of the detection method that, for Y = {r, A} and Y e = {r £ , A £ } (related to SM and 
LR respectively), reduces to 

| < Y £ > - < Y > | 
4c = , (12) 



y/l \a*{Y e ) + a*{Y)] 



where o-(Y £ ) 2 and <r(Y) 2 refer to the variances of Y £ and Y respectively (more details are in llllll ). 
Moreover we take into account the time needed to observe an escape by normalizing the index 
dxc w.r.t. the square root of the average ET (computed in the case of e = 0). Thus we have a 
rough measure of the amount of information per unit time extracted by the detection strategy. As 
clearly shown in Fig. @] in the most favorable condition, i.e. a signal frequency close to resonant 
frequency of the potential well, the index improvement for LR is as much as four-times that of the 
SM. The improvements are over the whole spectrum of the external signal, and in particular the 
performances increase at a frequency above the matching Eq.dlOl). while the two methods have the 
same efficiency below f2 = 0.2. With the LR the stochastic resonant peak of Eq.dTOl) disappears 
as expected [11]. The loss of information in the best frequency condition (for a low dissipation 
FP cavity) with respect to the matched filter is roughly a factor 2 at low Finesse values (J 7 = 2) 
and increases to a factor 3 at high values (J 7 = 1000). Putting it another way, the analysis of ETs 
requires about 4 9 times more data than the optimal matched filter in order to obtain the same 
detection performances. The number of escapes is an essential aspect of the detection of small 
signals (e « AV). When the signal is weak the number increases and thus constitutes a severe 
limitation to the detection. We have numerically checked that, at a reasonably high level of the 
index (fl2b . dxc — 8> the number of escapes increases as the square of the ratio between the noise 
intensity and the signal amplitude e, as for the washboard potential llllll . For instance to decrease 
the signal amplitude from e = 0.05 (as in Fig© to e = 0.0005 requires to increase the number 
of escapes from N ~ 400 to ~ 3 x 10 6 for the lower value of the Finesse (J 7 = 2), and from 
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N ~ 900 to N ~ 7 x 10 6 for the higher value of the Finesse (J 7 = 1000). 



0.14 
0.12 
O.l 

A 

^ 0.08 

V 

~p 0.06 
0.04 
0.02 



Driving Frequency O 

FIG. 4: Normalized Kumar-Carroll index as a function of the angular 
velocity of the applied signal for the SM approach (dashed line) and 
the Likelihood Ratio (solid line) keeping the same barrier height AV = 
0.07 for two Finesse parameters : T = 2 (< r >= 9.4, II m = 2.1 
and ui r = 2.9) and T = 1000 (< r >= 8.5, U M = 12 and u r = 6.6). 
The parameters of the system are: e = 0.05, <t>o = 0, D = 0.0125, 
7 = 10~ 6 . As a reference, the normalized cLkc index for the adapted 
filter is independent of the frequency and reads e/(2y/~D) — 0.22. 




VI. CONCLUSIONS AND OUTLOOKS 

We have shown that escape times can be effectively employed to determine the noise intensity 
and the presence of a coherent signal affecting a metastable system. To do so, we have focused on 
a FP interferometer, and we have explicitly discussed the advantages of the detection of the escape 
times: being mechanically decoupled, such characterization is in fact promising for measurements 
in the semiclassical or quantum regime. In spite of these advantages, had the escape sequence 
entailed a great loss of information respect to the full trajectory, the ET reduction would not be 
appealing. For this reason we have accurately compared the performances of the ET analysis with 
the full trajectory analysis (analysis performed through a matched filter, that can be demonstrated 
to be the optimal method). We have found that: 1) The noise intensity estimate is not distorted 
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(within 5%); 2) The noise estimate exhibits the optimal behavior N 1 / 2 for large N; 3) The esti- 
mates are more accurate in the limit of rare ETs; 4) The detection of a sinusoidal signal with the 
analysis of the average escape rate entails a loss of information (measured by the index dxc) of 
about an order of magnitude; 5) The detection of a coherent signal with a more accurate analysis 
based on a likelihood ratio estimator reduces the loss of the index dxc to a factor 2 4- 3; 6) Signal 
enhancement at a matching finite temperature given by Eq. (fT0l) . or stochastic resonance, is ap- 
parent in the analysis of the average escape time only. The more refined analysis with likelihood 
estimator enhances the performances that monotonically improve moving towards the resonant 
frequency, see FigjH 7) The detection sinusoidal signals exhibits the optimal behavior N 1 / 2 for 
large N, however the detection of small signals can be very challenging for mechanical systems, 
because the system is inherently slow and it might require an unpractical number of periods - 
as opposite to electronic systems such as Josephson junctions that are much faster and therefore 
potentially more suitable 113411 . 

These findings are of interest per se for the applications of Fabry-Perot pendular cavities in 
fundamental physics. However, a natural question that arises is about their generality for other 
systems that can be characterized by ETs. In other words, one could ask the following questions: 
does the loss of information (as measured by the variance of the estimate) entailed in the usage 
of ET depends upon the particular shape of the metastable potential? We lack a general answer, 
however one could notice that these features are also shared by the washboard potential associated 
to Josephson junctions lllin . and therefore we conjecture that they might be a general feature of 
estimates obtained through ETs. 

Moreover, there are two findings that are specific for ETs from FP potential, eq.©. First, 
the analytic evaluation of the average ET, eq.©, together with the assumption of an exponential 
distribution, captures the main features of the noise estimate with ETs, although with a discrepancy 
(around 30%) of the standard deviation, see Fig. |3] This finding is therefore useful to evaluate 
the performances of the estimates with ETs as a function of the FP parameter without heavy 
stochastic numerical simulations. Second, the nonlinearity parameter (the Finesse J 7 ) deteriorates 
the performances of FP as coherent signal detector, see Fig. HI thus giving a qualitative suggestion 
on FP design. 

We underline that in this paper we only illustrate the principle of operation of the ET detection 
through FP based interferometers, and we have shown the potentiality of the full analysis of the 
ET statistics. Essentially we suggest that the ET measurements method could be interesting if a 
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favorable balance between operational simplicity and information loss can be achieved. 
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